Electron spin decoherence in isotope-enriched silicon 
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Silicon is promising for spin-based quantum computation because nuclear spins, a source of mag- 
netic noise, may be eliminated through isotopic enrichment. Long spin decoherence times T2 have 
been measured in isotope-enriched silicon but come far short of the T2 = 2Ti limit. The effect of 
nuclear spins on T2 is well established. However, the effect of background electron spins from ever 
present residual phosphorus impurities in silicon can also produce significant decoherence. We study 
spin decoherence decay as a function of donor concentration, 29 Si concentration, and temperature 
using cluster expansion techniques specifically adapted to the problem of a sparse dipolarly cou- 
pled electron spin bath. Our results agree with the existing experimental spin echo data in Si:P 
and establish the importance of background dopants as the ultimate decoherence mechanism in 
isotope-enriched silicon. 
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Long electron spin decoherence times in silicon are of 
significant interest in producing low-error rates for quan- 
tum computation. Very long spin echo decay times T 2 
have been reported [l|-l3jj in isotope-enriched silicon (i.e., 
reduced nuclear spin concentration). Resource require- 
ments for quantum error correction are significantly re- 
duced as the qubit fidelity improves, which motivates bet- 
ter understanding of the limits of T 2 even with isotope 
enrichment. The ultimate decoherence time is theoreti- 
cally limited by inelastic decay mechanisms (spin-lattice 
relaxation) on a time scale of T\ . The increase of T 2 upon 
reducing the 29 Si nuclear spin concentration P, Q is now 
well understood 0, However, even the highest pu- 
rity Si wafers contain traces of dopant impurities, usually 
phosphorus, at levels ~ 10 12 — 10 14 cm~ 3 . Their electron 
spins are coupled by dipolar interactions, causing fluctu- 
ations that induce qubit spin dephasing. In this Letter, 
we develop the necessary theory to examine decoherence 
of a central spin in a sparse bath of nuclear and electron 
spins. We find excellent agreement with existing Si spin 
echo data showing that existing spin decoherence mea- 
surements in Si may already be limited by the coupling 
of the donor electron spin to the P donor spin bath rather 
than the Si nuclear spin bath. As a consequence, further 
isotopic enrichment, an extremely expensive procedure, 
may not provide any more advantage in the eventual con- 
struction of a Si spin quantum computer. In fact, we find 
that in the presence of donor-induced spin decoherence, 
T2 may actually increase when some 29 Si is present. 

We study here the central spin decoherence problem 
of a donor electron spin among spins of other donors 
and 29 Si. Because of coupling among the spins, a par- 
ticular donor electron spin will experience fluctuations of 
its energy splitting in a phenomenon known as spectral 
diffusion (SD). 29 Si- induced SD calculated using a clus- 
ter expansion technique [H, well approximated at the 



lowest order with independent contributions from each 
pair fj, is in excellent agreement with experiments for 
Si:P 0, i, 0j and Si:Bi 8] donors. With a firm founda- 
tion rooted in a precise quantum mechanical formulation, 
this was a significant advance over the long history of 
phcnomenological, stochastic models [9 11]. These pre- 
vious techniques [H, IB], however, are applicable to rela- 
tively dense and weakly-coupled spin baths and cannot 
accurately treat SD due to randomly located donors in 
which the strength of interaction to the central spin is 
no different than between bath spins; neither can they 
handle very low concentrations of 29 Si rigorously. A dis- 
joint cluster approach was applied to the relatively sparse 
bath of carbon spins for the SD of nitrogen- vacancy de- 
fects in diamond [12j [ . Exact numerics were applied 13[ in 
the central spin decoherence problem of dilute dipolarly- 
coupled spins. Our approach in this Letter is based upon 
the cluster correlation expansion (CCE) 14] that refor- 
mulates the cluster expansion technique 4j such that a 
large bath approximation is not necessary, making new 
regimes of the SD problem accessible. 

We consider an ensemble of Si:P donor electron spins 
over varied donor concentrations, Ce (for electron), and 
29 Si concentrations, Cjv (for nuclear). We use parts per 
million (ppm) of lattice sites for Cn ■ We include dipo- 
lar and hyperfine interactions among spins. Assuming 
a large applied magnetic field (100 mT is sufficient) in 
the z direction, we use an effective Hamiltonian in which 
Zeeman energies are conserved among the electron and 
nuclear spins independently, allowing only flip-flop dy- 
namics: rl = HE + H n + He-n, where 



He = 5>£ d ( R * - + § r§t - ^i^j], (1) 
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H N = ^ d ( r ™ - r rn)[%I- + - 4/*ft], (2) 
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with the dipolar interaction strength given by d(r) = [1 
3(r z /r) 2 }/4r 3 , and 



29 Si concentration, C N (cm 3 ) 



hi{K) = y|*i(R)| 2 



(3) 
(4) 
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written in atomic units; a factor of fi./(47reo) is implied for 
the Hamiltonian. The hyperfine interaction, hi(r), may 
be approximated by the dipolar interaction, d(r) when 
r is far outside the wave function of donor i. Electron 
spin operators are written as S with i or j indices and 
Ri position vectors. Nuclear spin operators are written 
as I with n or m indices and r„ position vectors. The 
gyromagnetic ratios of the electron and 29 Si nuclear spins 
are 7E = 1.76 x 10 n (T s)" 1 and lN = 5.31 x 10 7 (T s)" 1 
respectively. The wave function of each donor electron, 
\ti(r), is the Kohn-Luttinger wave function of a phospho- 



rus donor impurity in silicon, as described in Ref. 11|. 
In addition to the Hamiltonian-governed free evolution, 
we model spin echo refocusing pulses as ideal spin flips. 

To compute the decoherence time of a qubit in Si:P 
system, we take one of our donor electrons to be the "cen- 
tral" spin, say i = 0, and simulate a Hahn spin echo on 
that donor electron to remove the effects of static noise. 
Our dominant decoherence is due to flip-flopping bath 
spins: 29 Si-induced and donor-induced SD. We display 
agreement with experiment for over 5 orders of mag- 
nitude in Cat, Fig. [TJ maintaining agreement into very 
sparse densities. This Letter presents procedures we have 
developed to accomplish this substantial (and very com- 
putationally demanding) task. 



We previously |4j computed the decoherence for C 
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1000 ppm using a cluster expansion technique which 
works well for dense spin baths. For sparse baths, we use 
the CCE [3], applicable to both small and large spin 
baths, with some adaptations. The CCE has a simple 
and self-evident formulation which we now describe. We 
define C(t) = Pt±(t)/pti(0), the off-diagonal element of 
the reduced density matrix of our central spin after per- 
forming a spin echo sequence over the duration t — 2r, a 
refocusing pulse occurring at time r. The spin echo fig- 
ure of merit is the modulus of C(t). Next, for a given set 
(cluster) of electron or nuclear bath spins, C, we define 
Lc(t) to be the resulting C[t) when we only include flip- 
flop terms in our Hamiltonian [Eqs. <(T][3j)] that involve 
elements of C; all S^S* interactions are included in our 
implementation. Then, we recursively define 



L c {t) = Lc{t)l J] L c ,(t). 



(5) 
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By tautology, £(t) = Y\ c Lc(t), providing a way to break 
the problem into independent factors coming from each 
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FIG. 1: Decay times (T2 for the Hahn echo) of Si:P donor 
electron spins for various Cjv- At high Cjv, contact hyperfine 
interactions dominate and T2 is dependent upon the magnetic 
field direction relative to the lattice orientation. At low Cn, 
T2 is dependent upon Ce, and eventually dominated only 
by dipolar interactions (which includes dipolar-approximated 
electronuclear interactions). Experimental results are shown 
as square symbols, from Ref.Q, and a star symbol, from Ref.0. 



set of bath spins. At short times, the smallest nontrivial 
clusters dominate the decay; successively larger clusters 
become significant with increasing evolution time. 

These cluster expansions work well by perturbative ar- 
guments in the regime where the interactions among the 
bath spins are weak relative to the interaction with the 
central spin. Thus, 29 Si-induced SD is well-approximated 
when including only 2-clusters. Donor-induced SD, how- 
ever, is much more challenging because the interaction 
strengths among the bath spins and with the central spin 
are comparable. In fact, we find that if we compute the 
CCE expansion for different spatial configurations and 
different initial spin states, the average over these config- 
urations and states can diverge rapidly when we include 
4-clusters. We attribute this to the fact that different 
configurations of a sparse bath, or even different states of 
a given spatial configuration, can have very different con- 
vergence time scales for CCE. However, we find that the 
CCE is well behaved if we average over spin states within 
the CCE definitions, that is, C(t) = (p^(t)) j/(p^(0)) j, 
where J represents each spin state and p^ is calculated 
with J as the initial bath state. 

Of course, averaging over all spin states exactly would 
be prohibitively difficult. We find, however, that it is 
sufficient to average over spin states in the following self- 
consistent manner. Choose a spatial configuration and 
a spin state \ J) = ® n |j n ) that serves as a template for 
spin state variants. Let T be a set of clusters (e.g., up to a 
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FIG. 2: Spin echo resulting from decoherence induced by var- 
ious concentrations of (a) 29 Si (Cat), (b) background phospho- 
rus donors (Ce), or (c) the combination of both. In (b), static 
29 Si-induced Overhauser field variations between donors sup- 
presses their decoherence-inducing flip-flops, (a)-(c) use £ as 
a scaling parameter for both concentration and inverse time, 
(d) T2, defined as the 1/e decay time, for various Cn and 
Ce- The star symbol indicates an experimentally obtained 0] 
result at Cn ~ 50 ppm and is in good agreement with the 
corresponding 50 ppm theoretical curve. 



certain size) that we include to approximate the solution: 

^=n^ wc,r) . (6) 

eer 

where K(J,C,T) is the set of all spin states that may 
differ from J only for spins in superclusters of C that are 
contained in T. That is, 

fC(J,C,T) = {J'| 3C G T,C D C,T>(\J), \J')) C C'}, (7) 

where T>(\J),\J')) is the set of spins whose state differs 
between | J) and | J'). Then we define 



Lq - ( L c)Ke!c/ Yl 



(8) 
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where L J C solves the Lq problem for the given spin state 
J. Importantly, this yields the exact spin state average 
solution for £p in the limit that T includes all clusters 
(J becomes irrelevant). Furthermore, it may be com- 
puted relatively efficiently. With proper bookkeeping, 
each Hamiltonian (for a given cluster and external spin 
state) need only be diagonalized once, and each L J C need 
only be computed once and raised to the proper power 
to be multiplied into the solution. 

We use heuristics and cutoffs to determine the T set of 
clusters to include, trying to minimize the set necessary 
to approximate the solution well. We heuristically fa- 
vor clusters with strong interactions forming a connected 
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FIG. 3: (a)-(c) correlate with (a)-(c) of Fig. [2] respectively, 
but plotted in terms of error (one minus the spin echo) on a 
log-log scale in order to highlight the low-error behavior that 
is relevant for fault-tolerant quantum computation, (d) Plots 
TbflO -4 ], t ne time at which the echo error is 10 -4 , for various 
Cn and Ce- 



graph over the entire cluster and we employ cutoffs in the 
number of clusters, resonance energies, and distance from 
the central spin. We compute ensemble average results, 
such as shown in Figs. [TJ [2 and|3l by averaging results of 
different spatial configurations and J spin state templates 
for a given set of cutoffs. These cutoffs are adjusted until 
we obtain consistent, convergent results. 

We present, in Figs. [2] and [3j ensemble averaged spin 
echo results for varied Ce and Cat, both separately and 
combined, We use £ as a scaling parameter to illustrate a 
perfect correspondence between concentrations and in- 
verse time when decay is dominated entirely by 1/r 3 
dipolar interactions (Cat < 50 ppm). These results show 
behavior ranging from decay dominated by 29 Si- induced 
SD to decay dominated by donor-induced SD. When not 
dominated by 29 Si-induced SD, the presence of 29 Si can 
actually prolong coherence, see Figs.Q]and[2](d), because 
Overhauser field variations suppress donor flip-flops; a 
similar effect is noted [loj with respect to qubit concen- 
tration and layout. The initial decay, shown clearly in 
Fig. [31 behaves differently; in this regime, any beneficial 
effect of 29 Si is fairly insignificant. 

Computing or measuring ensemble averages, as shown 
in Figs. and 02 has limited utility in the scope of quan- 
tum computation. It is more informative to compute 
the full distribution of results that come out of the con- 
siderable sample-to-sample variation, which are signifi- 
cant especially for a central spin with dipolar coupling 
to a dilute bath [HI . Figure |4] addresses this for donor- 
induced spectral diffusion by showing error distribution 
information for each spin echo time independently; es- 
sentially, this gives a performance guarantee for various 
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Donor-induced spectral diffusion yields 
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FIG. 4: (top) Maximum error from donor-induced (without 
29 Si) SD for at least 50%, 90%, or 99% of the central spins due 
to random spatial variations of the bath spins, (bottom) Cor- 
responding fits of low-error behavior of 1 — exp [— (t/T q ) nq ] « 
(t/T q ) nq in the 10" 4 error regime. 



fractions of possible donors. At short times, the ensemble 
average echo decay error is actually dominated by statis- 
tical outliers as the top panel of Fig. 0] demonstrates. 
Used in the bottom panel of Fig. SJ we introduce T q and 
n q as figures of merit that characterize initial decoher- 
ence at short times, as appropriate for quantum infor- 
mation considerations. These are obtained by fitting the 
error to 1 - exp [-(t/T q ) n «] « {t/T q ) n » in the 1CT 4 er- 
ror regime (motivated by common fault-tolerance thresh- 
olds). These results have direct implications for quantum 
computer architecture designs and error analysis (16J. 

Apart from reducing Cg, donor-induced SD may be 
suppressed by polarizing the background donors ther- 
mally at temperatures that are readily achieved in spe- 
cialized refrigerators [17]. In the low-error limit where 
donor- induced SD often dominates [Fig. [3], the er- 
ror is proportional to the number of contributing 2- 
clusters which must have opposite spin polarization 
in order to flip-flop; thus, T q oc {p^Pi)~ nq , where 
p t/i = exp(±E z /k B T)/2cosh(E z /k B T) from Boltz- 
mann statistics with Ez as the electron Zeeman energy 
splitting corresponding to about 1.3 K per Tesla. 

To conclude, we adapted the cluster correlation ex- 
pansion [3], by retaining all Ising-like interactions and 
interlacing spin state averaging in a self-consistent man- 
ner, to study decoherence induced by a background of dy- 
namical donor electron spins in silicon. We demonstrate 
that approaching the T2 = 2T\ limit through isotopic en- 
richment in Si is impossible in the presence of a finite 
concentration of unpolarized donors. Unavoidable donor 
impurities in the background make this limit, where T\ 
of one hour has been reported at 1.25 K 18j, unattain- 
able, though prospects improve if the electrons may be 



thermally polarized. While the presence of some 29 Si can 
actually increase T2 considerably by suppressing donor- 
induced decoherence, this effect is fairly insignificant in 
the short time (low-error) regime important for quantum 
computation. We introduce T q and n q to describe deco- 
herence at short times and discuss the effect of statistical 
variation of impurity locations on decoherence. Variation 
in the decoherence of different donors becomes extremely 
significant in the regime of low impurity concentration, 
a crucial consideration for designing quantum computer 
architectures and determining fabrication requirements. 
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